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ABSTRACT 

Many engineering structures, such as airplane wings, beams and shafts are subjected to higher torsional forces 
today due to advancement in Structural Engineering, in terms of size and technology. In this paper, we analyzed the 
resistance of circular beams, of different engineering materials, to their corresponding twisting moments. We obtained the 
torsional rigidity for the different beams as the ratio of twisting moment to the angle of twist per unit length. It is observed 
that torsional rigidity of the beams is a function of their areas and the engineering material they are made up of. 
Specifically it is observed that the circular beam made up of brass engineering material has the greatest torsional rigidity 
among the twelve engineering materials considered. 
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INTRODUCTION 

When a beam is transversely loaded in such a manner that the resultant force passes through the longitudinal shear 
central axis, the beam only bends and no torsion will occur. When the resultant force acts away from the shear central axis, 
then the beam will not only bend but also twist. [1, 2, 4] 

Torsion is twisting about an axis produced by the action of two opposing couples acting in parallel planes [5]. 
Another name for couples is torque or twisting moment. Torsional rigidity of a beam is a ratio of moment to the angle of twist 
per unit length [6]. When torsion is applied to a structural member, its cross-section may warp in addition to twisting. If the 
member is allowed to warp freely, then the applied torque is resisted entirely by torsional shear stresses (called St. Venant' s 
torsional shear stress). If the member is not allowed to warp freely, the applied torque is resisted by St. Venant' s torsional 
shear stress and warping tension. This behavior is called non-uniform torsion [1, 2, 3, 4]. 

Beams of non-circular section tends to behave non-symmetrically when under torque and plane sections do not 
plane. Also the distribution of stress in a section is not necessarily linear[12]. 

St. Venant' s theory is usually applied when the cross-section is non-deformable out of its plain or those deformations 
are very small [10]. 

Consider a circular beam with length / , with one of its bases fixed in the xy — plane, while the other base (in the 

plane Z = / ) is acted upon by a couple whose moment lies along x — axis. The beam twists through an angle determined by 
the magnitude of applied couple and the modulus of rigidity of the beam. The amount of twist produced can thus be used to 
determine the applied force [5]. 
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Saint- Venant (1885) was the first to produce the correct solution to the problem of torsion of bars subjected to 
moment couples at the ends 

In material science, shear modulus or modulus of rigidity, denoted by jU or G , is defined as the ratio of shear stress 
to the shear strain [6]. 

FORMULATION OF THE PROBLEM 

For a beam of constant circular cross-section subjected to torsion, the St. Venant's torsion is given by [1, 2, 3, and 4] 

T sv = Imf- (1) 

dz 

Where, 

<p is the angle of twist (twist angle), 
jU isthemodulus of rigidity, 

Tgy is St. Venanttorsion, 

/ is the polar moment of inertia, 

Z isthedirectionalong axis of themember. 

df 

= Twist rate 

dz 

By symmetry, anysection of the beam perpendicular to z — axis remains perpendicular to this axis during 
deformation and the action of the couple will merelyrota teeach section through some angle (f) , called the angle of twist [5]. 

Theamount of rotation wil lclearlydependonthedistance of the section from the base Z = 0, and since the deformations are 
small, the amount of rotation 0 is proportional to the distance of the section from the fixed base. 

Theangle of twist can then be written as 

f = az, (2) 
Where a isthe twist per unitlength. 

It is therelative angular displacement of a pair of crosssectionsthat are unitdistanceapart. Letw be the 
displacementalongz - axis. 

w = 0, if the crosssection of the beam remainplaneafter deformation. Since w is independent of z , we write 

w = af\x,y) (3) 

Where ^isthetorsion function 
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Consider a particle originalty at (x,y,z) , since/ = az. , for the displacement of this particle 

u = - azy , v — azx , w = 0. 

The angle of twist / can also be written as [5] 



f-r- 



Where^T is the applied torque and / is the length of the beam. 




Figure 1: Twisting of Circular Section 

i 2 




P ° (x v x 2 ) 
P'° (r, + u u x 2 + u 2 ) 
Figure 2: Showing the Displacement Along x x and x 2 

Assumptions 

• The bar is straight and of uniform cross section. 

• The material of the bar has uniform properties. 

• The only loading is the applied torque which is applied normal to the axis of the bar. 

• The bar is stressed within its elastic limit. 



(4) 
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The Eulerian strain is given in terms of displacement u , v , w by [5] 
1 



c = —ill. . + u . ) 

u 2 hJ h 



Where u, , - — '- , i = 1, 2, 3; / = 1, 2, 3 



u x = u,u 2 = vandw 3 = w 



x l = x, x 2 = y andx 3 = z 



P Uy — u — - azy — - ax 2 x- i — -fx 2 



(5) 



u 2 = v = - azx = ax^^ = fx Y 

The representative strains, in matrix form, are as follows [5] 



0 
0 



0 
0 



f -{- x 2 + ^i) 



2 
/ 



0 



Where 

fx, " fx 2 
The stress-strain relation is given as [5] 



i E , u 

d - = ( e a + : — — e kAj) 



1^2 

fx, 



1 + u 



1 - 2u 



Where dL isthe stress tensor, £ isthe Young' smodulus and M isthePoisson's ratio. 
Sub stituting the respectivestrains in the stress-strainrelation, weget 

? o o tfC#V * 2 )S 

0 0 Jt/C^ + Xj): 

gxa(T;i - Xj) xf(f*,2 + x x ) 0 ^ 



(6) 



(7) 



(8) 



Where X is a lameconstant, given as 
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E 

x = (9) 

2(1 + u) 

Recall that a function of X j and X^ ■ 

\ for/ = l,^n = 0 

fori = 2,^22 = 0 

P f^U + ^22 = 0(10) 

This is the Laplace equation. Hence, ^is a harmonicfunction. 

Since^/^is a harmonicfunction in thesimplyconnected region R, representingthecross-section of thebeam, there 
exist sananalyticfunction^+ iy of the complex variable X + iy where)? (x,y) is a harmonic conjugate of 
ThefunctionssatisfytheCauchy-Riemann equations, namely 

- 11 • (i.e.,^1 = y a ) (11) 



fx 1 fx 2 



And 



f 2 ^ f 2 y 



(i.e., f 1 ! = - y,) (12) 



fx 2 fx, 
The harmonic conjugate y of 3^is given as 

y = -(x 2 + x\) on cf 
and y (13) 

N 2 y = Oin R 

This is a Dirichlet problem. 
ANALYSIS 

The torsional rigidity of a beam is defined as a ratio of moment to the angle of twist per unit length [5]. The torsional 
rigidity of a beam with circular cross-section is given as [3, 5] 

D = -y 1 (14) 
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WhereM 3 is the resulting momento on the surfacex 3 = / and is given as [8]: 

M 3 = qq (x - x 2 d 3l )dA (15) 

R 

where J 32 an d cL^ are stress tensors with j = xf g X 2 ~i"and = x f 6 1" 

i ( Hx 2 0 l^lxj 0 

On the surface X 3 = / , we havex^ = (x^,X 2 , 0) , hencex 3 = 0 . Therefore, M j = 0 and M 2 = 0 

P M 3 = M = 0C)( x i4s2~ x 2 d3i)dA (16) 

R 

Letusconsidertheharmonicfunction [8] 

y = c 2 {x\ - x 2 ) + k 2 (17) 
where C , k are constants. 

0 0 O 0 1 O 0 

f> C (x{ - X 2 ) + k = ~( x i + x 2 ) ontheboundary, or 

(-- c 2 )x 2 + (-+ c 2 )x 9 2 = & 2 (18) 
2 2 

The curve defined by this equation is an ellipse 



2 2 
Xi X 9 

■7+ -f = 1 (19) 

.2 l2 



a 2 ^ 



2 1 tC • k 

If we choose c < — and a = — , , £? = . _ , then 

2 JH 



2 la 2 - Z? 2 2 aV 
c = , k = 



2 a 2 + b 2 ' a 2 + b 2 



la 2 - b 2 , , aV 



2 a 2 + b 2 ' 1 a 2 + b 2 

- 2xfa 2 x 7 
So, flu , becomes — 

a z + a z 
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2xfb 2 x, 

Similarly, a 32 becomes L 

a 2 + b 2 

From equation (16); the torsional momen to becomes 



2xf c - 
M = 2 2 P 2 QO JC i 2 ^ + a 2 QQxldA± (20) 



R R 



2x1 (a 2 I x + b 2 I x ) (21) 



a 2 + b 2 

where/ and / are the moments of inertia of the elliptical section about thex, - and X 0 - axes. 

pab 2 pa 3 b 

But/ = and / = , so we have 

x i 4 x 2 4 

pxfa 3 b 3 

M = — (22) 

a 2 + b 2 

For a beam with circular cross-section; 

a - b - radius (r) 

= pxfa 3 a 3 = pxfa 4 = pxfr 4 
a 2 + a 2 2 2 

The torsional rigidity of the circular beam can be written as 
„ M pxr 4 

D = — = (24) 

/ 2 

butA - pr 

P> — = r 2 (25) 
A 2 4 

-7 = r 4 (26) 
P 

Substiruting equation (26) into (24), wehave 
pxA 2 

D = (27) 

2p 2 
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^ xA 2 

D = (Torsional rigidity) (28) 

2p 

P J — = X = (29) 

A 2 2(1 + u) 

A 2 F 

V D = A n \ (30) 
4p(l- v) 

For a circular beam; the moment of inertia about thex - axis is given as 
3 

l * = (31) 

Also, the moment of inertia about the y - axis isgiven as 
3 

P r 

K = *V (32) 

Thus, the polar moment of inertia about thex 3 - axis is given as 

3 3 

pr pr 

7 p = V 4 (33) 

3 

pr 

l P = ~Y ( 34 ) 
Hence, St. Venanttorsion becomes 
p_r\ndq_ 

1 SV T J ' 

2 dx 3 
With the boundary conditions 
/(x 3 ) = 0atx 3 = 0,/(x 3 ) = / at x 3 = /; 

To find the twist angle J , we integratea long the length of thebeam as shownbelow 

6f^3 = 6 (36) 

ax 3 ml 

p 

IfT S(/ , Ifl and I p A — nr 2 constants along the beam. 
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ml p 



\ f - ^SV X 3 

ml p 



If length of the beamalongz - axis is/ , then 



T I 
_ 1 sv L 

ml „ 



From equation (3), the displacement is given as 
x 3 = af\x,y) = ai\x 1 ,x 2 ) 
Displacement along X j and X 2 axes are 



Uy — ~ ClX yK 2 and U 2 — (XX ^ X ^ 



respectively , since displacement along^ 3 - iszero. 
If the length of the beam alongXj - axis is/ 
Mj = - f X 2 &w\U 2 = f%\ (since/ = CLZ ) 



U y = 



Tz Tz 
ml p ml p 



If the length of the beam along X 3 - axis is/ , then 
77 77 



Mj = 



ml 



ml , 



(37) 



(38) 



(39) 



(40) 



(41) 



(42) 



(43) 



From equation (13), the harmonic conjugate y , of the torsion function which is a function of X ^ and X 2 , is 
considered for the different values of X j and X 2 and plotted on a graph as shown in Figure 2. 

Table 1: The Following Chart Gives Typical Values for the Modulus of Rigidity, Young' S 
Modulus and Poisson Ratios for Different Engineering Materials [6, 9] 



Engineering 
Materials 


Modulus of 

Rigidity 
(Psi X 10 6 ) 


Young's 
Modulus 
(Psi X 10 6 ) 


Poisson 
Ratio 

(w) 


Beryllium copper 


6.7 


17 


0.285 


Brass 


5.8 


102-125 


0.331 
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Table 1: Contd., 


Bronze 


6.5 


96 - 120 


0.34 


Copper 


6.58 


17 


0.355 


Iron (Malleable) 


9.4 


28.5 


0.271 


Magnesium 


2.39 


6.4 


0.35 


Molybdenum 


17.16 


40 


0.307 


Monel 


9.57 


26 


0.315 


Nickel silver 


5.6 


18.5 


0.322 


Nickel steel 


10.8 


29 


0.291 


Titanium 


5.94 


27 


0.32 


Zinc 


6.1 


12 


0.331 



For this paper, we consider a circular beam, for twelve different engineering materials of diameter 1 .2m, angle of 
twist of 30°, length of 10m and the Torque (T) as the St. Venant Torsion (T sv ). Table 2 shows the calculated values of polar 

moment (/ ), St. Venant Torsion (T ' S y ) and torsional rigidity ( D ). 

i.er = 0.6m ,/ = 30°, / = 10m , T = T sv 

Table 2: Calculated Values of Polar Moment/ , St. Venant Torsion T sv and Torsional Rigidity D 



Engineering 
Materials 


m 


p 


T 

1 sv 


D 


Beryllium copper 


6700000 


0.34 


6834000 


1.347 


Brass 


5800000 


0.34 


5916000 


8.683 


Bronze 


6500000 


0.34 


6630000 


8.207 


Copper 


6580000 


0.34 


6711600 


1.278 


Iron 


9400000 


0.34 


9588000 


2.283 


Magnesium 


2390000 


0.34 


2437800 


0.483 


Molybdenum 


17160000 


0.34 


17503200 


3.116 


Monel 


9570000 


0.34 


9761400 


2.013 


Nickel silver 


5600000 


0.34 


5712000 


1.425 


Nickel steels 


10800000 


0.34 


11016000 


2.287 


Titanium 


5940000 


0.34 


6058800 


4.757 


Zinc 


6100000 


0.34 


622000 


0.918 



Displacement along x { and x 2 for circular beam of different engineering materials, with/ = 10m , 

I p = 0.34 

Table 3: Displacement Along x j and x 2 for Circular Beam with Different Engineering Materials 







u x 


u 2 


1 


i 


-30 


30 


2 


2 


-60 


60 


3 


3 


-90 


90 


4 


4 


-120 


120 


5 


4 


-150 


150 


6 


6 


-180 


-180 


7 


7 


-210 


210 


8 


8 


-240 


240 


9 


9 


-270 


270 


10 


10 


-300 


300 
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The Relationship between jU , T sv and D 

When two variables X and y are related, they are said to be correlated [1 1]. In order to define the linear relationship 

and the amount of linear relationship between jU , T sv and D , we adopt the concept of regression and correlation 
coefficient. We deduced the following: 

ju = 2970474 + 0.688458r sv (44) 

and 

R = 0.878477 (45) 
Similarly, we have 

D = -4.9xl0~ 9 + 3. 104177// (46) 

And 

R = -0.00659 (47) 
Where r is the correlation coefficient, ( —1 < r < 1). 
Recall that 

D = Torsional rigidity, 

/ = Polar moment of inertial about I, - axis, 
P 3 

T sv = St. Venant torsion, 

tn = Modulus of rigidity 

The Relationship between the Cross-Sectional Area (A) and Torsional Rigidity (D) 

In this section we are interested in finding out if there is any relationship between the cross sectional area of the 
beams and torsional rigidity of the beams of different engineering materials. As shown in tables (4)-(8), where r is the radius , 
we considered five cases where the diameter of the cross sectional area is given values 1.2, 2.2, 3.2, 4.2 and 5.2 respectively. 

Table 4: Torsional Rigidity of Beams with Different Materials When Cross-Sectionalareais 1.130973m 2 





r 


A 


E 


V 


D 


BerylliumCopper 


0.6 


1.130973 


17 


0.285 


1.346606406*10 A 6 


Brass 


0.6 


1.130973 


113.5 


0.331 


8.679859374*10 A 6 


Bronze 


0.6 


1.130973 


108 


0.34 


8.203776869*10 A 6 


Copper 


0.6 


1.130973 


17 


0.355 


1.277040024*10 A 6 


Iron 


0.6 


1.130973 


28.5 


0.271 


2.282412788*10 A 6 


Magnessium 


0.6 


1.130973 


6.4 


0.35 


4.825486315*10 A 5 


Molybdenum 


0.6 


1.130973 


40 


0.307 


3.115152316*10 A 6 


Monel 


0.6 


1.130973 


26 


0.315 


2.012530533*10 A 6 


Nickel Silver 


0.6 


1.130973 


18.5 


0.322 


1.426568663*10 A 6 


Nickel Steels 


0.6 


1.130973 


29 


0.291 


2.286475953*10 A 6 
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Table 4: Contd., 


Titanium 


0.6 


1.130973 


27 


0.32 


2.082019130*10 A 6 


Zinc 


0.6 


1.130973 


12 


0.331 


9.176943822*10 A 5 



Table 5: Torsional Rigidity of Beams with Different Materials When Cross-Sectionalareais 3.801327m 2 



Engineering Material 


R 


A 


E 


V 


D 


BerylliumCopper 




3.801327 


17 


0.285 


1.521270401*10 A 7 


Brass 




3.801327 


113.5 


0.331 


9.805696062*10 A 7 


Bronze 




3.801327 


108 


0.34 


9.267862436*10 A 7 


Copper 




3.801327 


17 


0.355 


1.442680787*10 A 7 


Iron 




3.801327 


28.5 


0.271 


2.578457226*10 A 7 


Magnessium 




3.801327 


6.4 


0.35 


5.451384652*10 A 6 


Molybdenum 




3.801327 


40 


0.307 


3.519208725*10 A 7 


Monel 




3.801327 


26 


0.315 


2.273569408*10 A 7 


Nickel Silver 




3.801327 


18.5 


0.322 


1.609166175*10 A 7 


Nickel Steels 




3.801327 


29 


0.291 


2.583047410*10 A 7 


Titanium 




3.801327 


27 


0.32 


2.352071 148*10 A 7 


Zinc 




3.801327 


12 


0.331 


1.036725575*10 A 7 



Table 6: Torsional Rigidity of Beams with Different Materials When Cross-Sectionalareais 8.042477m 2 



Engineering 
Material 


R 


A 


E 


V 


D 


BerylliumCopper 


1.6 


8.042477 


17 


0.285 


6.809505981*10 A 7 


Brass 


1.6 


8.042477 


113.5 


0.331 


4.389222715*10 A 8 


Bronze 


1.6 


8.042477 


108 


0.34 


4.148477786*10 A 8 


Copper 


1.6 


8.042477 


17 


0.355 


6.457723385*10 A 7 


Iron 


1.6 


8.042477 


28.5 


0.271 


1.154168245*10 A 8 


Magnessium 


1.6 


8.042477 


6.4 


0.35 


2.440147154*10 A 7 


Molybdenum 


1.6 


8.042477 


40 


0.307 


1.575267148*10 A 8 


Monel 


1.6 


8.042477 


26 


0.315 


1.017694453*10 A 8 


Nickel Silver 


1.6 


8.042477 


18.5 


0.322 


7.202944776*10 A 7 


Nickel Steels 


1.6 


8.042477 


29 


0.291 


1.1 56222903 *10 A 8 


Titanium 


1.6 


8.042477 


27 


0.32 


1.052833377*10 A 8 


Zinc 


1.6 


8.042477 


12 


0.331 


4.640587891*10 A 7 



Table 7: Torsional Rigidity of Beams with Different Materials When Cross-Sectional Area is 13.85442m 2 



Engineering 
Material 


R 


A 


E 


V 


D 


BerylliumCopper 


2.1 


13.85442 


17 


0.285 


2.020751239*10 A 8 


Brass 


2.1 


13.85442 


113.5 


0.331 


1.302521396*10 A 9 


Bronze 


2.1 


13.85442 


108 


0.34 


1.231 079267 *10 A 9 


Copper 


2.1 


13.85442 


17 


0.355 


1.916358186*10 A 8 


Iron 


2.1 


13.85442 


28.5 


0.271 


3. 425045690* 10 A 8 


Magnessium 


2.1 


13.85442 


6.4 


0.35 


7.241245398*10 A 7 


Molybdenum 


2.1 


13.85442 


40 


0.307 


4.674675445 *10 A 8 


Monel 


2.1 


13.85442 


26 


0.315 


3.020053630*10 A 8 


Nickel Silver 


2.1 


13.85442 


18.5 


0.322 


2.137505955*10 A 8 


Nickel Steels 


2.1 


13.85442 


29 


0.291 


3.431 142979*10 A 8 


Titanium 


2.1 


13.85442 


27 


0.32 


3.124329957*10 A 8 


Zinc 


2.1 


13.85442 


12 


0.331 


1.377115133*10 A 8 
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Table 8: Torsional Rigidity of Beams with Different Materials When Cross-Sectional Area is 21.23717m 2 



Engineering 
Material 


R 


A 


E 


V 


D 


BerylliumCopper 


2.6 


21.23717 


17 


0.285 


4.748200693*10 A 8 


Brass 


2.6 


21.23717 


113.5 


0.331 


3.060561278*10 A 9 


Bronze 


2.6 


21.23717 


108 


0.34 


2.892692239*10 A 9 


Copper 


2.6 


21.23/ 1 / 


1 H 
1 / 


0.355 


a cftiftn/;iftft*iftAO 
4.502906190*10 o 


Iron 


2.6 


21.23717 


28.5 


0.271 


8.047900209*10 A 8 


Magnessium 


2.6 


21.23717 


6.4 


0.35 


1.701490304*10 A 8 


Molybdenum 


2.6 


21.23717 


40 


0.307 


1.098418091*10 A 9 


Monel 


2.6 


21.23717 


26 


0.315 


7.096282046*10 A 8 


Nickel Silver 


2.6 


21.23717 


18.5 


0.322 


5.022541645*10 A 8 


Nickel Steels 


2.6 


21.23717 


29 


0.291 


8.062227127*10 A 8 


Titanium 


2.6 


21.23717 


27 


0.32 


7.341302272*10 A 8 


Zinc 


2.6 


21.23717 


12 


0.331 


3.235835711*10 A 8 



RESULT DISCUSSIONS 

The numerical calculations were carried out for a circular beam of length 1, with one of its bases fixed in the xy- 
plane, while the other base (in the plane z= 1) is acted upon by a couple whose moment lies along the z-axis (x 3 -axis). As an 
illustration, the length 1 of the beam is taken to be 10m, the diameter, 1.2m, the angle of twist, 30° and the polar, 0.34. (i.e 
l=10m, r=0.6m, ^=30°, I p = 0.34. Twelve different engineering materials with different values of E, raand v were considered. 
The results are shown on the various tables and figures. It is observed from table 2, that circular beams of brass engineering 
material has the highest torsional rigidity under St. Venanttorsion, while circular beam made of magnesium engineering 
material has the lowest torsional rigidity. Clearly from the value of R in equation (45), which is approximately 0.9, it shows 
there is a strong positive correlation between the modulus of rigidity and St. Venant torsion, and equation (44) shows there is a 
linear relationship between modulus of rigidity and St. Venant torsion. Also, that the higher the modulus of rigidity of the 
engineering material the higher the St. Venant torsion. However, the value of R in equation (46), which is approximately - 
0.01, shows that apart from the fact that torsional rigidity and modulus of rigidity have a negative correlation, the correlation 
is very weak. This implies that the value of the modulus of rigidity has a little negative effect on the resistance of twist of a 
circular beam. The effect of the cross sectional area of the circular beam is shown in tables (4),(5),(6),(7) and (8). It can easily 
be seen that the cross-sectional area of the beams affect the torsional rigidity of the beams. The wider the cross-sectional area 
the higher the torsional rigidity. 

CONCLUSIONS 

This work deals with the analysis of torsional rigidity of circular beams with different engineering materials 
subjected to St. Venant torsion. The torsional rigidity of these beams were calculated as a ratio of twisting moment to the 
angle of twist per unit length. It is shown that the circular beam made of brass engineering material has high torsional rigidity, 
relatively, when subjected to St. Venant' s torsion. Concerning the cross-sectional area of circular beams made of different 
materials; it was deduced that the wider the cross sectional area the higher the torsional rigidity and vice versa. This 
phenomenon is of great importance, especially in the field of civil and mechanical engineering. 
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